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ABSTRACT 


The potential for remote sensing of ocean wave direction us- 
ing passive polarimetric microwave observations was investigated. 
A fixed-beam 91.65 GHz polarimetric radiometer was mounted 
on the NASA DC-8 aircraft during the Tropical Ocean Global 
Atmosphere/Coupled Ocean Atmosphere Response Experiment 
(TOGA/COARE January - February 1993). Several experiments 
were performed during which the DC-8 was flown in constant 
bank-angle turns at ^ 1.5 km altitude to obtain azimuthal scans 
of the sea surface at fixed observation angles Data at 65 from 
nadir are consistent with previous findings using 19- and 37-GHz 
Special Sensor Microwave Imager (SSM/I) satellite observations 
and support the claim that a broadband emission mechanism is 
responsible for the azimuthal brightness signatures. Accordingly, 
a tihed-facet geometrical optics (GO) model of the surface vs as 
developed to investigate emission from deterministic and random 
striated surfaces. Laboratory measurements of polarimetric emis- 
sion at 92 GHz from small-amplitude water-waves corroborate 
this model. 


INTRODUCTION 


Recently there has been much interest in passive remote sens- 
ing of ocean surface wind and wave direction. An investigation 
using satellite data from the SSM/1 has shown that the surface 
wind vector can be determined using dual-polarized brightness 
temperatures (Wentz, 1992). Observations of the upwelling ra- 
diation over wind-driven oceaD surfaces showed brightness tem- 
peratures that were related to wind direction. In this paper we 
describe a model for these azimuthal brightness signatures and 
discuss the potential of using fully polarimetric radiometry for 
ocean wind direction sensing 

It is proposed that a broadband mechanism is responsible 
for the upwelling brightness temperature signatures reported here 
(91.65 GHz) and from the 19- and 37*GHz SSM/1 data (Wentz, 
1992). Tbe brightness signatures can be broken down into contri- 
butions from the first and second azimuthal Fpurier harmonics. 
We explain the presence of each Fourier harmonic as due to ocean 
wave asymmetry (the first harmonic) and geometrical scattering 
(the second harmonic). 

A tilted-facet geometrical optics (GO) model of the surface 
emission incorporating geometrical optics is used to explain tbe 
physical origin of the second harmonic in the azimuthal bright- 
ness signatures. Laboratory measurements of a sinusoidally stri- 
ated surface corroborate the GO model and demonstrate its appli- 
cability in predicting polarized- thermal emission from relatively 
smooth water surfaces. A Monte- Carlo analysis of simulated one- 
dimensionally rough ocean surfaces using the GO model correctly 
predicts the phase of the second azimuthal harmonic in observa- 
tions over the open ocean. Although the firsl-barmonic is not 
modelled in this analysis it is hypothesized to be due to ocean 
wave asymmetry and foam. 


TOGA/COARE OBSERVATIONS 


/d 



To characterize the polarimetric upwelling azimuthal bright- 
ness temperature at 92 GHz over an open ocean several exper- 
iments were performed during TOGA/COARE in January and 
February 1993. A fixed-beam polarimetric 92 GHz radiometer 
capable of measuring 7V, 7* and Tv * 2 Rc(E v El) was installed 
on tbe NASA DC-8 airborne platform. Severs! constant bank- 
angle turns were performed at ~ 1.5 km altitude to obta n 3o0 
azimuthal scans of the sea surface at fixed observation angles. 
Nadir video imagery and aircraft navigational information was 
mied to orovide sea surface and wind truth. 


An azimuthal scan at an observation angle 6 t = 65® (Fig 1) 
showed variations in 7V of ^1.3 K and 7V of ^2.5 a 

nounced upw'ind/downwind asymmetry in T v and T h (Kunkee and 
Gasiewski, 1994). Data from this scan are consistent with previ- 
ous findings using 19- and 37-GHz SSM/1 satellite data (Wentz. 
1992), suggesting that the upwelling azimuthal signatures are the 
result of broadband emission. 



Figure 1: Measured azimuthal variation of 7V and TV (91.65 GHz) 
with wind direction observed at $ $ — 65® over open ocean under 
light wind conditions during TOGA/COARE. 
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The above brightness temperature signatures can be conve- 
niently described using the azimuthal Fourier amplitudes, B on 
and phases, $ on . defined by: 

T on {6,) = £ B m (9.) cos[n«J, + *.„(«.)] - m6.) (1) 

«=0 


where T (c) is the calm water brightness, a is the polarization, 
n is the^azimuthal harmonic and is the azimuthal observa- 
tion angle with respect to the water wave vector (wave azimuthal 
angle). For an upwind observation 4>, = 0°. Both first-order 
(n = 1) and second-order (n = 2) variations can be found in the 
TOGA/COARE data at 6 , = 65°. For the vertical brightness, 
T ^ 0.29cos (<t>, 4 69°) 4 0.13 cos[2(^ - 36°)] and for the hori- 
zontal brightness, 7V ~ 0.81 cos(d>, 4 47°) 4 0.53cos[2(<£, - 42°)]. 
Note that the brightness temperature minimum is not directly 
in the downwind direction with an average phase shift of -48 
shown. This shift is due to primarily to the presence of multiple 
wave trains and light winds. 

The first-order azimuthal harmonic (£ a i) is hypothesized to 
result from wave asymmetry and increased amounts of oce^i foam 
on the leeward side of the wave. Leeward-side foam will be par- 
tially hidden from the view of the radiometer during downwind 
observations resulting in a lower observed brightness temperature 
compared to the upwind observation. The first-order harmonic 
mav also be due to ocean wave asymmetry: the leeward side of 
the wave is generally steeper and will thus appear brighter due 
to reflection of background radiation incident on the surface from 
the warmer near-grazing angles. 


GEOMETRICAL OPTICS MODEL 


The second-order harmonic from a striated water surface 
can be explained using geometrical optics (GO). A GO model o 
the thermal emission from a sinusoidally striated water surface 
was developed to investigate the azimuthal variation in the up- 
welling polarimetric brightness temperature (Gasiewski and Kun- 
kee, 1993). Using the GO model the surface is composed of spec- 
ularly reflecting facets, each contributing to the total upwelhng 
radiation. The model incorporates multiple scattering up to third 
order. 

Use of the GO model is justified at 92 GHz by considering 
the domain of applicability of the Kircboff approximation (KA). 
Because KA is often the most direct and least time consurmng 
method for predicting surface emission, it is unfortunate that 
there seems to be no unanimous decision regarding its applica- 
tion to surfaces of varying roughness. Wirgin (1983), claims that 
the KA is applicable to a sinusoidal surface at normal incidence 
for h/A < 0.011(A/A), where h and A are the surface amplitude 
and wavelength respectively, and A is the electrical wavelength. 
Perhaps a more universal criteria involves the minimum local ra- 
dius of curvature of the surface (Beckmann and Spizzichino 1963), 
this region of validity is min(2tr c cos 0) > 1 where r c is the local 
radius of curvature of the surface, k is the electrical wave num- 
ber and 0 is the local angle of incidence. Thus, the GO model is 
expected to degrade at observation angles near-grazing, however, 
the above criteria are easily met at 92 GHz over a wide range of 
water surface roughness and observation angles. 

Indeed, based on laboratory measurements of a sinusoidally 
striated water surface, Wirgin 1 * criteria appears conservative for 
studies of polarimetric emission from water surfaces. To show 
this, direct measurements of the first three upwelling Stokes pa- 
rameters (TV, TV and Tv) over a rotatable water wave tank were 
performed at several observation angles using a precision cali- 
brated polarimetric 92 GHz radiometer (Gasiewski and Kunkee, 
1993). The measurements were used to corroborate the GO model 
and to establish relationships between T v , TV and Ty and wave an- 
gle. All experiments were performed under clear- skv background 
conditions; coincident radiosonde measurements of atmospheric 
temperature and water vapor profiles were obtained for comput- 
ing sky background brightness temperatures. 


Comparisons of the GO model with measurements are shown 
in Figs. 2 and 3 for a sinusoidally striated surface described by 
z = h sin(2?rx/A) with h/A = 0.05 and 0, = 65°. The re- 
sults verify the applicability of GO with multiple scattering for 
calculating polarimetric thermal radiation from relatively smooth 
curfarpc at near-Brewster angle observation. 



Azimuthol Wave Angle (degrees) 


Figure 2: Measured and calculated T v and TV brightnesses over a 
sinusoidal water surface for 0 , = 65°, T w = 288A , h/A = 0,05 
and SRH = 50%. 



Azimuthol Wove Angle (degrees) 


Figure 3: Measured and calculated Tv brightnesses over a sinu- 
soidal water surface for 0. = 65°, T w = 288 A', h/A = 0.05 and 
SRH = 50%. 

For a sinusoidal surface the GO model shows that the second- 
harmonic Tv brightness variations are in phase quadrature with 
the T v and Th second-harmonic magnitudes (Figs. 2, 3). Due 
to this phase relationship, coincident observations of Tv with 
and 7V will facilitate the retrieval of wave and ocean surface wind 
direction by reducing a four-quadrant direction ambiguity to a 
two-quadrant direction ambiguity. 


MONTE-CARLO OCEAN WAVE ANALYSIS 


A wind-driven sea surface is more complex than a sinusoidal 
model and must be described statistically. For the open ocean, 
the slope distribution has variance (Cox and Munk, 1954; Wilheit, 

1979 ) : 

el(w) = 0.003 4 0.0048tr (2) 

where w is the wind speed in m/s at a height of 20-m above the 
surface. To incorporate ocean wave statistics, a Monte- Carlo anab 
ysis was performed using a set of one-dimensional wave-covered 
surfaces and the GO model. The simulated ocean surfaces have 
wavenumber spectra with a power law of slope -3.0 and a slope 
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variance c\(l m/s) = 0.037. Results for 6, = 20° and frequency 
f - 12 GHz (Figs. 4, 5) show azimuthal brightness variations of 
magnitude and phase close to those obtained by Johnson using a 
diffraction model (Johnson. 1993). 



Figure 4: Calculated TV and 7* brightnesses frorn a simulated 
ocean surface for 6, = 20' , T u = 300 K, as 0 037. 



Figure 5: Calculated T v brightnesses from a simulated ocean sur 
face for 6, = 20' . T, = 300 K. o\ 0.037. 


Results at 92 GHz asuming 507c surface relative humidity 
(SRH) and clear-air conditions are shown in Fig. 6 for an obser- 
vation angle of 65°. From these data * o2 = J' 80 "^ 7 ) “ d . pp 
This is consistent with as computed from the 1 UUA/tUAKt 
data (after translation by a - -48° phase correction) (Fig. 1). 
The calculated values of B q2 are several times larger than the 
TOGA/COARE observations. However, this model is somewhat 
incomplete in that ocean foam, two-dimensional slope spectra and 
cusping are not yet considered. 



DISCUSSION AND CONCLUSIONS 

Passive polarimetric observations at 92 GHz of the sea sur- 
face under light sea conditions show variations of -1.3 h and -2.5 
K in the vertically and horizontally polarized brightness temper- 
atures over an azimuthal .can of 360” at a constant observation 
anele of 65°. An upwind/downwind brightness asymmetry in T v 
and T h was observed in this data supporting the claim that oceaD 
surface wind direction can be unambiguously observed using pas- 
sive polarimetric remote sensing. 

Unambiguous recovery of surface wind/wave direction on the 
open ocean will rely on the first-order azimuthal terms defined 
in (1) As expected for sinusoidal surfaces (symmetric surfaces) 
all fi are zero for r odd n, however, since real ocean waves are 
not symmetric this does not hold in general. The combination 
„{ B , 4 0 (an upwind/downwind asymmetry) and T V id phase 
quadrature with T v and T k will provide a way in principle to 
eliminate all ambiguities in the azimuthal angle of the retrieved 
o cean surface wind/wave direction. 

The amplitude and phase of the azimuthal variations are 
consistent with satellite observations at 19 and 37 GHz from 
the SSM/I. This consistency suggests that a broadband emission 
mechanism is responsible for the ocean wave azimuthal signatures 

The geometrical optics model is shown to be valid in applica- 
tions to water wave remote sensing at 92 GHz. Not only does the 
region of validity extend to near-Brewster observation angles, but 
results that are consistent with a diffraction model are uso ob- 
tained at lower frequencies using a statistical ocean surface. Th 
GO model also correctly predicts the phase of the second_order 
harmonic in the TOGA/COARE (92 GHz) observations. Better 
agreement with the observed second-order harmonic magnitude 
and accurate estimates of the first-order harmonic magnitude and 
phase are expected using a more complete statistical ocean surface 
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